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The finite duration of tlie collisions in Fermionic systems as expressed by the retardation time in 
non-Markovian Levinson-type kinetic equations is discussed in the quasiclassical limit. We separate 
individual contributions included in the memory effect resulting in (i) off-shell tails of the Wigner 
distribution, (ii) renormalization of scattering rates and (iii) of the single-particle energy, (iv) colli- 
sion delay and (v) related non-local corrections to the scattering integral. In this way we transform 
the Levinson equation into the Landau-Silin equation extended by the non-local corrections known 
from the theory of dense gases. The derived nonlocal kinetic equation unifies the Landau theory of 
quasiparticle transport with the classical kinetic theory of dense gases. The space-time symmetry 
is discussed versus particle-hole symmetry and a solution is proposed which transforms these two 
exclusive pictures into each other. 



I. INTRODUCTION 

The generalization of the Boltzmann equation (BE) 
towards dense interacting quantum systems is a still de- 
manding and unsolved task. In dense systems, the mean 
time between successive collisions becomes comparable to 
the collision duration. In this case all particles cannot be 
described as a system of weakly interacting quasiparticles 
but a set of particles bound in two - particle correlations 
has to be accounted for in a distinct manner. In the spirit 
of chemical reactions, one can introduce the quasiparticle 
density, (a density of single-atom molecules), and the 
density of particles in the correlated states, ric (twice the 
density of bi-atomic molecules). The total density n is 
their sum 

n = Uf -\- ric- (1) 

In parallel with the partial pressure of molecules, V = 
kBT{nf -\- ^Uc), and the Guldberg-Waage law of mass 
action law, ric = ^"-/j even a small fraction of the corre- 
lated density, ric nf, contributes to the second virial 

coefficient, V ~ ksT (^n — ^Kn'j^. The correlated den- 
sity thus plays an important role in the thermodynamic 
behavior of the system. 

While the equation of state of the ideal gas, V = ksTn, 
does not reflect any microscopic properties of particles 
in the system, the second virial coefficient directly fol- 
lows from the particle-particle interaction. In the early 
time of statistical mechanics, the experimental second 
virial coefficient has been used to deduce interaction po- 
tentials. Surprisingly, the firmly established concept of 
the equilibrium virial expansion has never found a corre- 
sponding position in the theory of non-equilibrium sys- 
tems although a number of attempts have been made to 
modify the BE so that its equilibrium limit will cover 
at least the second virial coefficient The achieved 



corrections to the BE have a form of gradient or nonlo- 
cal contributions to the scattering integral. For a hard- 
sphere gas the non-local correction is trivial, therefore the 
main theoretical focus was on the statistical correlations 
P-|l5| . It turned out that the treatment of higher order 
contributions is far from trivial as the dynamical statisti- 
cal correlations result in divergences that are cured only 
after a re-summation of an infinite set of contributions. 
Naturally, this re-summation leads to non-analytic den- 
sity corrections to the scattering integral |ll|-|l5|. The 
moral of these hard-sphere studies is that beyond the 
non-local corrections one must sum up an infinite set of 
contributions; the plain expansion leads to incorrect re- 
sults. From this point of view, the approximate statistical 
virial corrections implemented by , although possibly 
reasonable, are not sufficiently justified by the theory. 

Real particles do not interact like hard spheres, in par- 
ticular when their do Broglie wave lengths are compara- 
ble with the potential range. An effort to describe the 
virial corrections for more realistic systems resulted in 
various generalizations of Enskog's equation . 
By closer inspection one finds that all tractable quantum 
theories deal exclusively with the non-local corrections. 
The statistical correlations in quantum systems would 
require an adequate solution of three-particle collisions 
(say from Fadeev equations) which are now intensively 
being studied |^. A systematic incorporation of three- 
particle collisions into the kinetic equation, however, is 
not yet fully understood, therefore we discuss only binary 
processes. 

An alternative way to describe the correlation has been 
developed for Fermi systems at very low temperatures. 
Since the Pauli principle excludes all but the zero-angle 
scattering channels, the correlation can be re-cast into a 
renormalization of the single-particle energies known as 
the quasiparticle energies. While quasiparticles behave 
in many aspects like free particles, the density depen- 
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dence of the quasiparticle energies results in non-trivial 
thermodynamic properties. In accordance with the fo- 
cus either on zero-angle or finite-angle scattering, in the 
literature one can find two distinct classes of quantum 
kinetic equations. 

The zero-angle class are equations determining the 
time evolution of the momentum and space dependent 
quasiparticle distribution function f{k,r,t). This leads 
to kinetic equations of the Landau-Silin type 



df dedf dedj < > 



(2) 



Here the drift term is characterized by the quasiparti- 
cle energies e{k,r,t) and the collision integral consists of 
scattering-in and -out terms which are summarized in the 
selfenergy as functionals of the quasiparticle distribution 
function, (T^[/]. The collision rates are reduced by the 
wave function renormalization, which is given by the fre- 
quency derivative of the real part of the selfenergy a at 
the pole value z^^ — ^~ ^ \uj=e- This kinetic equation is 
local in time and space, i.e., it describes the collisions as 
instantaneous point-like events. Similarly to the BE, the 
local collisions yield no virial corrections. The virial cor- 
rections are thus covered exclusively by the quasiparticle 
energy while no correlated density appears. 

The finite-angle class are equations which include the 
finite duration of collisions, i.e., they have non-Markovian 
scattering integrals. These kinetic equations are usu- 
ally developed for the reduced density matrix or the 
Wigner function r, t). The Wigner distribution obeys 
a Levinson-type of kinetic equation [^ , ^ 

dpt de^^ dpt de^^ dpt 
dt dk dr dr dk 



2Im / dtGi 



(l-Pt)S? 



< 



(3) 



Here the drift term is essentially determined by the 
Hartree-Fock mean field, = k"^ /2m + a^f, and the 
double-time selfenergies S^-^ are given as functionals of 
the Wigner distribution. Explicit time arguments of the 
collision integral, t < t, show the retardation which re- 
quires to include the propagation, G^^, from the retarded 
time t to the time t of balancing changes. Having the 
non-Markovian form, this kinetic equation is capable to 
describe finite duration effects and corresponding virial 
corrections. 

The Landau-Silin and Levinson equations have com- 
mon features with the original BE. In particular, they 
express the balance between the drift given by the gradi- 
ent terms on the left hand side and the dissipation given 
by the scattering integral on the right hand side. On the 
other hand, there is a remarkable difference in the actual 
physical contents of these common ingredients, namely 



P ¥^ /> £ 7^ ^-nd the right hand sides of (||) and ^ 
are rather different. 

In the present paper we argue that the differences 
between the phenomenological quasiparticle picture, 
Eq. (|), and the BBGKY-type]] equation (|) can be char- 
acterized in terms of the retardation of the collision in- 
tegral of (^). We will show that there are various contri- 
butions to the total retardation, each responsible for dif- 
ferent features seen in the quasiparticle picture extended 
by non-local corrections. 

If both approaches are equivalent, the Levinson equa- 
tion (^) must contain terms on the collisional side which 
should be possible to rearrange into gradients such that 
the quasiparticle energies of the drift side in the Landau- 
Silin equation (||) are obtained. The mutual relation be- 
tween the Wigner distribution p and the quasiparticle 
distribution / is of essential importance for such a rear- 
rangement. 

In Sec. II we show that a part of the retardation of 
the Levinson equation (||) describes the off-shell motion 
of particles. This off-shell motion can be eliminated from 
the kinetic equation which requires to introduce an effec- 
tive distribution (the quasiparticle distribution /) from 
which the Wigner distribution p can be constructed Mq) 



P = .f + 



duj 



P 



d 



2tt uj — e duj 



{{i-f)<y^-fK 



(4) 



This relation is the extended quasiparticle picture derived 
for small scattering rates. The limit of small scattering 
rates has been first introduced by Craig Ml . An inverse 
functional f[p\ has been constructed in |45|1 . For equi- 
librium non-ideal plasmas this approximation has been 
employed by [^ , ^ and under the name of the gener- 
alized Beth-Uhlenbeck approach has been used by ji^ j 
in nuclear matter for studies of the correlated density. 
The authors in have used this approximation with 
the name extended quasiparticle approximation for the 
study of the mean removal energy and high-momenta 
tails of Wigner's distribution. The non-equilibrium form 
has been derived finally as the modified Kadanoff and 
Baym ansatz [Q. We will call it extended quasiparticle 
approximation. 

We show that the retardation is responsible for gra- 
dient terms by which the mean-field drift of the Levin- 
son equation (^) differs from the quasiparticle drift of 
the Landau-Silin equation (j^). In Sec. Ill we discuss 
the remaining parts of the retardation, a piece which is 
compensated by the decay of internal propagators, and a 
piece which describes the collision delay. Special empha- 
size is put on internal double counts of correlations in the 



'Equations derived from the Born-Bogoliubov-Green- 
Kirkwood-Yvon (BBGKY) hierarchy of reduced densities. 
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Levinson equation. Related to this question we present 
in appendix C a discussion of a common met pitfall in 
literature. In Sec. IV we complete the kinetic equation 
with the finite collision delay and provide related non- 
local corrections in space. Our main result is the nonlo- 
cal kinetic equation derived here in an alternative way to 

dSg dfa 

dr dk 



dfa _^ dsg dfa 
dt dk dr 



= E 



V< faKil - fa){l - fb) - fafbil - f'a){l - fi) 

(5) 



with - faik-q-AK,r~A3,t~At), fi - fbip + q- 
Ak, r-A4, t-At), fa = faik, r, t) and fb = fb{p. r-Aa, i). 
The differential cross section V ~ is proportional 

to the square of the amplitude of the T-matrix. All non- 
local corrections are given by derivatives of the scattering 
phase shift (j) — Iva. Int^iuj, k,p,q, t, r) Q , 



d<j) 1 90 ld(f) 



dk' 



Ao 



dp dq dk ' 



d(j) d(t) 



We discuss in Sec. IV the space-time symmetry of this 
collision integral and link the correlated density with the 
collision delay. In Sec. V we summarize and refer to nu- 
merical solutions of the derived nonlocal kinetic equation. 
The most important formulae for the T-matrix are pre- 
sented in appendix A and the technique how to derive 
nonlocal corrections from gradient expansion is summa- 
rized in appendix B. 



- i (Go 'G< - G<Go^) = - i (E«G< - G<E^) 

+ i (G^S< - S]<G^) . (7) 

The retarded and advanced functions are introduced as 
5^(1,2) = -ieih - t2){B> + B<) and B^{1,2) = 
i9{t2 — ti){B^ +B^). The products are understood as in- 
tegrations over intermediate variables (time and space). 
The Hartree-Fock drift term has the form 



Go 1(1, 2) 



dh 2m J ^ ' 



(8) 



The Kadanoff and Baym equation ^ is more general 
than any of the above kinetic equations which can be 
derived from it with the help of specific approximations. 
One possibility is to rearrange the terms such that they 
yield the Landau-Silin kinetic equation for the quasipar- 
ticle distribution function ( ^5| ) when the gradient expan- 
sion is applied |5^,|5^,Q. The other possibility is to col- 
lect the terms on the right hand side of (^ into a non- 
Markovian collision integral which results in the Levin- 
son equation (||) for the Wigner distribution [p7[-p9[. 
The connection between both equations, however, has 
remained dubious because in the Levinson equation the 
drift term is controlled exclusively by the Hartree-Fock 
selfenergy while in the Landau-Silin equation the drift 
term is represented by the full quasiparticle energy. In 
the following we demonstrate how the Levinson equation 
can be converted into the Landau-Silin equation. To this 
end we first derive the Levinson equation from (^. 



A. Precursor of the Levinson equation 



II. KINETIC EQUATIONS 

First we want to show that the Levinson equation can 
be transformed into the kinetic equation of Landau-Silin 
type and vice versa. Therefore we first derive the kinetic 
equation for the Wigner distribution from the real-time 
Green function technique. For a homogeneous system 
similar investigations have been presented by Bornath at 
al Q . Here we extend their approach to inhomogeneous 
systems and want to focus specially on the physical con- 
tents of both equations. This approach will furnish us 
with a set of Kramers-Kronig relations needed to link 
the Wigner and the quasiparticle distributions. 

We consider the two independent correlation func- 
tions for Fermionic operators G>(1,2) ~ (a(l)a+(2)) 
and G^(l,2) = (a"'"(2)a(l)), where cumulative variables 
mean time and space, 1 = i, x. The time diagonal 
part of G^ yields the Wigner distribution, p{xi,X2,t) = 
G^(l, 2)tj For the correlation function G^ the 

Kadanoff and Baym equation of motion reads Q 



On the time diagonal, ti,2 — t, the Kadanoff and Baym 
equation (0) yields an identity 

-t{Go'G< -G<Go'){t,x,,t,x2) 
t 

= j dt' Jdx'(G>{t,xi,t',x')i:<{t',x',t,X2) 

— oo 

+j:<{t,xi,t',x')G>{t',x',t,x2)^ 

t 

- J dt' Jdx'(G<{t,xi,t',x')J:>{t',x',t,X2) 

— OO 

+E> {t, XI , t', x')G< (t', x', t, X2)) . (9) 

On the left hand side there is the Hartree-Fock drift and 
the right hand side contains a non-Markovian collision 
integral. Note that correlations beyond the Hartree-Fock 
field are exclusively in the collision integral. 
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1. The Wigner representation 

Now we concentrate on the gradient expansion of the 
cohision integral. The quasi-classical approximation is 
achieved using the Wigner mixed representation for space 
arguments 

^ ( , X1+X2 ti+t2 ^ 

o- I K, ^ , — ^ — ,ti - <2 

d{xi -X2)e-*'=("i-"^)E(ti, a;i, i2, X2). (10) 

As a rule, we use the lower case to denote operators in 
the full Wigner representation (momentum, space, time, 
frequency) and the hat to denote the Wigner representa- 
tion in space and double-time representation. Note that 
time arguments of & now denote the center-of-mass time, 
, and the difference time, ^1—^2- These two Wigner 
representations are identical for p and a which do not 
depend on the difference time. 

In the representation (^o|) the identity (||) reads 



d k d \ 
at m or J 

+2sin Q(a,^5f - dld^)y'^\k,r,t)p{k,r,t) 

00 

= e^^(^-{dldi-dld'^)^ jdr 



/ T T 

.[g>{k.r.t~ -,T)a<{k,r,t- 
+&<{k,r,t - ^,T)g>{k,r,t - ^, -t) 

T T 

-g<{k,r,t- -,T)a>{k,r,t - 

T T 

-a>{k,r,t-- -,T)g<{k,r,t- -, -t] 



(11) 



Here, the superscripts 1 , 2 denote that the partial deriva- 
tives apply only to the first and second function in the 
product, e.g., {d^d^^fab=p^p^. 



2. Gradient approximation 

The power expansion of the goniometric functions in 
( prj ) defines the expansion in space gradients. This ex- 
pansion goes to infinite order, but in the following we 
will restrict our attention to the linear expansion. It is 
customary to abbreviate the linear gradient terms with 
the help of the Poisson brackets, 

„ . /l/.^r^!P ^1 r^9C\ , r ,T dadh da dh 
2sm [^-idldl - dld^;)y ^ {aM = -- - --. 

(12) 



The linear approximation in space gradients is straight- 
forward, however, one has to be careful about time argu- 
ments of functions in the collision integral of (jl^). For 
example, the first and second product seem to form an 
anticommutator in which the linear gradients in space 
cancel. This is not true because of time arguments. 

The gradient approximation in time requires a special 
treatment because of the lower integration limit. Due to 
this limit the integral does not define a standard matrix 
product with respect to the time variables. We will give 
this expansion an explicit notation. 

In equilibrium, the functions g< and d< do not depend 
on the center-of-mass time. For slowly evolving systems 
the center-of-mass dependence is smooth and weak. Ac- 
cordingly, in the collision integral of (|ll|), we can expand 
the center-of-mass time dependence around the time t in 
powers of r. Since all functions in (11) have the same 
center-of-mass time, t — -j, it is possible to write the lin- 
ear expansion of (^l|) with respect to time gradients in a 
compact form 



(13) 



with e^^ = ^ + a^^ . The zero order gradient term, 
1 = 1^ — /<, reads 



00 

I<^^l+'-(d^dl-dld^)^ Jdr 



X (<7> {t, T)g< {t, -r) + g< (<, r)a> {t, -r)) . (14) 

is obtained from via the interchange of particles 
and holes, >^<. In the first order gradient term, R ~ 
_R> — i?< , is given by 



X {a>{t,T)g<{t,-T)+g<{t,T)a>{t,-T)). (15) 

Now we are ready to turn all functions into the Wigner 
representation for time and energy. 



{k,r,t) = / dre"^^ a{k,r,t,T). 



(16) 



We will write the energy argument as the subscript and 
keep other arguments implicit in most of the formulas. 

After a substitution of cr's and g's of the Wigner rep- 
resentation (|l^) into (^4|), one can integrate out the dif- 
ference time T. The integration over time results in the 
6 function which describes processes on the energy shell, 
and the principle value terms which represent off-shell 
processes. 



y'dre'^"'-"'')^ =7r(5(cj- 



uj) + i- 



P 



(17) 
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The i5-function can be readily integrated out leaving both 
functions with identical energy arguments. The principle 
value part,p, is an odd function of energies u and w', 
therefore its non-gradient contributions cancel but the 
linear gradients survive. Formula (p^ thus turns into 



dujduj' 



P 



{a>,.g<}. (18) 



In the correlated part we express the time integral as 



d 



duj J 



= —inS'luj — uj') -] 



(19) 



The meaning of 5' and p' is seen from comparison with 
(p^). The S' and p' are odd and even functions in w — w', 
respectively. The gradient expansion of (^) thus reads 



dmduj' p' > < f dcj . > 



(20) 

The term R contributes to the kinetic equation ( p^ ) 
as a gradient correction linear in time. Since the second 
term of i?< in ( pO| ) is proportional to space gradients, we 
neglect this term leading to the second order contribution 
in gradients. 



B. Connection between the Wigner and 
quasiparticle distributions 

At very low temperatures, dissipative processes are 
known to vanish due to the Pauli exclusion principle. It 
is desirable to reorganize the kinetic equation so that the 
scattering integral will vanish in this limit too. As can 
be seen from (po|), the term dtR remains finite even for 
very low temperatures. It is possible to formally remove 
dtR if we shift this term on the drift side of ( |l3| ) and 
introduce a new distribution function. 



f = p-R. 



(21) 



We have denoted the new function as / to anticipate that 
it is indeed the quasiparticle distribution as will be shown 
now. 

Since we want to arrive at the kinetic equation for /, 
it is advantageous to write relation ( pT] ) in the opposite 
way so that we obtain p as a functional of f , p = f + R. 
Using ( pO| ) and the particle- hole conjugated term, R^ , 
one finds 



P = f- 



duj'duj p' 
(27r)2 uj' -u 



(22) 



This relation provides the so called extended quasiparti- 
cle approximation of the Wigner distribution. 



1. Wave-function renormalization 

To separate the on-shell and off-shell contributions of 
the correlation term, we use the imaginary part of the 
selfenergy, 7 = + and the spectral function, a — 
+ to rewrite ( p^ ) as 



/ 



dio'dw 



P 



(27r)2 oj' -oj 
With the help of the Kramcrs-Kronig relation 



(23) 



duo' 



P 



. I' 



duj 



2n UJ — LU 

where a is the real part of the selfenergy, we obtain 



(24) 



P = f + 



duj dau. 
2tt duj 



■9u + 



dio'dio 



P' 



(27r)2 uj' -UJ 



(25) 



The simplest quasiparticle approximation, ~ 2ttS{uj — 
e) and g^ — f 2TTd{uj — e), in the correlation terms yields 



, , da- 



/ + 



duj' 



p' 



2tt u' 



(26) 



In ( |2q ) the on-shell (quasiparticle) contribution is re- 
duced by the wave-function renormalization, 



z = l+ — 
doj 



(27) 



and the term with cr< provides the off-shell contributions. 

The exact renormalization is = 1 — d^^a, therefore 
( p7| ) is only the linear approximation in d^jcr. From the 
Kramers-Kronig relation ( p^ ) one can see that this corre- 
sponds to the linear approximation in the scattering rate 
7. We will keep all terms only up to this order in 7. 



2. Extended quasiparticle approximation 

From the relation between the Wigner and the quasi- 
particle distributions ( ^ ) one can deduce an approxi- 
mative construction of the correlation function g'^ as 
a functional of the quasiparticle distribution /. Since 
p = J j^g"^, we can write ( p6| ) as 



duj 



j '^[f z2.5{uj - e) 



UJ ■ 



(28) 



We have multiplied / with the 5 function and included 
this term into the integrand. 

One can see that the correlation function in the ex- 
tended quasiparticle approximation l5q,pQ], 



9^ = f z2tt5{uj - e) 



P 



UJ — £ 



(29) 
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satisfies (|2q). Perhaps it is not necessary to recall that 
the function / has been introduced here merely to sup- 
press the non-dissipative part of the collision integral, 
dtR- In |5^,|5^, the function / has been defined as the 
factor of the singularity of which justifies its inter- 
pretation in terms of the quasiparticle distribution. As 
relation ( p8| ) shows, both approaches lead to identical re- 
sults. 

In parallel with (|29|), for the hole correlation function 
one finds 



{I - f)z2TT5{Lo ~ e) + 



P' 



(30) 



From the spectral identity, a — + g"^ , it follows that 
( p9| ) and ( ^0|) are consistent with the extended quasipar- 
ticle approximation of the spectral function 



z 2-11 5{u} — e) 



P' 



uj — e 



(31) 



It is worthy to remark that (|31|) fulfills the spectral 
sum rule IBtI 



duj 
2^ 



a{k, UJ, r,t) = 1 



as well as the energy weighted sum rule |Q 
—uja{k,uj,r,t) = — +a^i{k,r,t). 



(32) 



(33) 



The latter one is violated in the simple quasiparticle pic- 
ture. 



C. Landau-Silin equation 

So far we have treated only the time gradients finding 
that the identity (13) can be expressed as 



dt 



(34) 



Now we rearrange this identity into the Landau-Silin 
equation for /. 

The time-local remainder / of the collision integral still 
includes the space gradients, / = — /y. The non- 
gradient part [the first term of (p^j is the Boltzmann- 
type scattering integral of the Landau-Silin equation 

Ib=J g(.9>'x< - g<a>) = . ((1 - /)a< - fa>) 

(35) 

In the second step we have used (^^. The off-shell con- 
tributions to (|3^) have not been neglected but cancel 
exactly. 

The gradient part which is the second term of ( p^ ) 
reads 



(27r)2 uj' -to 



dhJ 

W) 



(36) 



Here ( p^ ) has been used to remove the w'-integration for 
(T^ -|- (T^ and a similar Kramers-Kronig relation for the 
w-integration over the spectral function + .9^ has been 
applied. The function g = Heg^'^ ~ \{9^ + 9^) the 
off-shell part of the propagators. 

The mean-field drift {e^^,p} in (|3^ ) is not compatible 
with the quasiparticle distribution in the time derivative 
and the scattering integral. Moreover, the space gradi- 
ents from the collision integral, /y, have to be accounted 
for. Our aim is now to show that all gradients can be 
collected into the quasiparticle drift, 



{e^',p} + /y -{£,/}. 



(37) 



The quasiparticle energy differs from the mean-field en- 
ergy by the pole value of the selfenergy, e = e*^^ + (Tg . 

We will proceed in two steps. First we observe that the 
right hand side of (^^ includes only the on-shell contri- 
bution, therefore we show that the on-shell part of the 
left hand side equals the right hand side. Second, we 
show that the off-shell part of the left hand side of ( |37| ) 
vanishes. The on-shell parts of (29) used in ( ^6| ) result 
into 



= j dLo{a^,zf5{Lo-e)} 

= We, zf} - a'^ie, zf} + 2/{cr^, e} 

- {a„ zf} + (1 - z){e, zf} + zf{z ~ 1, e}. (38) 

Now we add the on-shell part of the commutator 
{e'^^,p}™ — {e^^, zf}. Using a rearrangement 

{e'^\zf} = {s~a,zf} 

= {£, f} + {z- l){e, /} + f{e, z} - {a, zf} (39) 
we obtain 

{e'^f , p}™ + = {£, /} - {£, (1 - zff}. (40) 

The last term is of higher order than linear in the damp- 
ing and can be neglected which confirms the relation ( |37| ) 
already from the on-shell parts. 

It remains to prove that all off-shell contributions to 
the left hand side of (|^), O = {e^\p]°^+I°« cancel, 
i.e., = 0. From the ofT-shell term of ( p9|) we find 

duj 



hf 



+ fT.,^^T<|-{5^,a<}j. (41) 



The term with e'^^ results from {e^^ , p}, while the others 
from /y . In the last term we use the extended quasipar- 
ticle approximation of the real part of the propagator 
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(42) 



which directly resuhs from ( |31| ) and the Kramers-Kronig 
relation. Using {:;^,o-<| = -{e, ^j^cr^} we can 
rewrite O as 




The linear expansion in the vicinity of the pole, <Tui—(J^ 
[ijj - e){z - 1) + 0(7^), yields 



O 



27r \ uj — e 



(44) 



The product (z — l)a^ is of second order in 7, i.e., the 
off-shell contribution O is negligible within the assumed 
accuracy. This completes the proof of relation (p7|). 

In summary, the requirement that the scattering inte- 
gral vanishes at very low temperatures directly leads to 
the concept of quasiparticles represented by relation ( |2^ ) 
between the Wigner and the quasiparticle distributions 



/ + 



duj 



P 



d 



271 Lu ~ e duj 



((1-/)^ 



(45) 



The space gradients of the collision integral renormalize 
the mean-field drift into the familiar quasiparticle drift 
( ^ ) . The resulting kinetic equation for the quasiparticle 
distribution has the structure of the Landau-Silin equa- 
tion (0) 



df de df de df 
dt dk dr dr dk 



z{{\~f)af~fa>). (46) 



We conclude this section that the Levinson type of 
equation (H) is equivalent to the Landau-Silin equation 
(^ up to second order in the damping 7 or in the ex- 
tended quasiparticle picture. 



III. ADDITIONAL RETARDATION OF THE 
LEVINSON EQUATION 

The link between the Levinson and Landau-Silin equa- 
tions shows that the collision integral of the Levinson 
equation cannot be fully interpreted in terms of the scat- 
tering processes but it includes various renormalization 
features. In the above treatment we have constructed 
a closed equation for the quasiparticle distribution func- 
tion while the reduced density is given by an additional 
functional of this quasiparticle distribution. There arises 
an important question whether it is not possible to con- 
struct the correlation functions, g^'^, directly from the 



Wigner distribution so that identity ( |11| ) turns into a 
closed equation for the Wigner distribution. To show 
the problems and drawbacks with this opposite way we 
discuss this construction for a homogeneous system. 

An approximation giving the Levinson equation is the 
GKB ansatz which for the homogeneous system reads 
[r>0] 

9^[k,t~ ^,t) = ig^{k,t - ^,T)p{k,t-T), 
g<{k,t-^, -r) = -ig^ik, t~^, -r)p{k, t-r). (47) 

For g^ one substitutes 1 — p for p. Note that this approx- 
imation includes the additional retardation between the 
center-of-mass time ^ — § and the time argument of the 
Wigner distribution which is t — r. This additional re- 
tardation causes many problems ranging from instability 
of numerical treatments over double counts of renormal- 
izations or virial corrections to incorrect interpretations 
of various correlation phenomena. Here we show how to 
handle this retardation within the linear approximation. 

If one assumes that the propagator depends only on 
the difference time r, i.e., it does not depend on any time 
dependent external fields or mean fields, it is possible 
to follow the approach developed above with only minor 
modifications. To make our discussion specific we limit 
our attention to homogeneous systems and employ the 
Galitskii-Fcynman T-matrix approximation of the self- 
energy. 



(7<(fc,t 







■9 (p,t 



-r) 



(27r)6- 

j dTiidTA'T^{k,p, q, Ti^)T^{k,p, q, ~ta) 



X g< [ k - q,t - -{th + TA + t),t - tr + TA 



X g< (^p + q,t-^(TR + TA+ t),t -TB. + taJ ■ (48) 

This selfenergy describes the process in which two par- 
ticles of initial momenta k ~ q and p + q are scattered 
into final states k and p. The function T^'^ is the re- 
tarded/advanced scattering T-matrix, see appendix A. 
Due to a finite duration of the scattering process, the T- 
matrices bring yet additional retardation times tr^a- We 
will show that retardation times tr^a are responsible for 
the virial corrections to the density. 



A. Decay of propagators and double counts 

First we focus on double counts following from the 
additional retardation resulting from the GKB ansatz. 
Since these problems appear already for weakly coupled 
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systems, it is sufRcient to restrict our attention to the 
Born approximation, T^'^ = V5{tb,,a)- The Levinson 
equation then reads 



d_ 

dt 



p(fc,t) = 2Re j ^Vlj drSit-^,T)D{t-r). 



(49) 



All distributions are collected in 

D{t) = {l-p{k,t)){l-p{p,t))p{k~q,t)p[p + q,t) 

- p(k, t)pip, t){l~p{k - t)){l-p{p + q, t)), (50) 

all propagators are collected in 

5(t,T) =5^(fc,t,r)5^(p,i,T) 

X g^{k-q,t,-T)g'^{p + q,t,-T). (51) 

We have used that the Wigner distribution is a real 
function and the complex conjugacy of propagators, 
g^{k, t, t) = [g^{k, t, — t)) , to reduce the expressions. 

Let us assume for the moment that the propagators do 
not depend on the center-of-mass time t. In this case the 
only t dependence in the collision integral of the Levinson 
equation (|9|) is due to distributions. Since for the Born 
approximation all distributions are retarded in an equal 
way, we can follow a modification of the approach used 
above. 

For slow perturbations one can expand ( p9| ) to the low- 
est order in the memory, D{t — t) = D{t) — TdtD{t), 
so that the collision integral splits into the Boltzmann- 
like scattering integral and the gradient correction, dtp — 
I + dtU, with 



^^'^W2Re|dr5(r), 



oo 



(52) 



(53) 



Please note that TZ and X are different from (^3|). Again, 
from (49) we can define a new function to remove the 
retardation. 



J^ = p-TZ, 



(54) 



and try to redefine the transport in terms of J-. This will 
us lead now to unavoidable double counts. 



that the de-phasing is the dominant part leaving the 
decay aside. Within the quasiparticle approximation, 
g^{k,T) w — ie""^'-'^, which is free of the decay, the time 
integrals over propagators yield 



oo 

2 Re y drSqp (r) 2tt5 {sk + Sp - Sk- 



OO 

2Rey" drr5'qp(T) = -2 



^^ ^k — q ^p-\-q 



(55) 



(56) 



After approximation (p5|), I from (|5^) turns into the 
Boltzmann-type scattering integral Ib of the Landau- 
Silin equation. Naturally, the on-shell contributions are 
identical since they do not depend on the retardation. 
Using (p3) in (53) one finds 



11 = 2 



dpdq 



VlD{t) 



P' 



£k—q £p+q 



(57) 



The function TZ is twice the off-shell contribution R to 
the Wigner distribution. To show this, we express (|5^) in 
terms of the Born selfenergy ( ^ ) and the approximative 
correlation functions, g^ = 2TrS{Lu — e)p, as 



7^= -2 



dw'duj 



(27r)2 LU-LU^ 



7 (^^'5^ - '^^'9u) • (58) 



Comparing with the second term of ( |22| ) one can see that 
the additional retardation results in the double count of 
the correlated part, TZ = 2R. Due to this double count, 
the function JF = p — TZ = p — 2R = f — R cannot 
be interpreted as the distribution of excitations because 
it has large regions of negative values. Please note a 
common met pitfall in literature related to this double 
count given in appendix C. 

We note that in numerical treatments of the Levinson 
equation neglecting or underestimating the decay, leads 
to a numerical instability of the equation, as it has been 
reported by Haug [ pO| . Other numerical solutions show a 
continuous increase of kinetic energy, e.g. figure 4 of |6l| ] 
or pulsation modes |6^] . The instability of the solution is 
a problem which shows that the Levinson equation is a 
difficult tool to handle. The double count of the off-shell 
contribution resulting in the negative effective occupa- 
tion factors is likely one of the reasons of this problem. 
A detailed discussion and numerical comparison of the 
Levinson equation with the Kadanoff and Baym equa- 
tion can be found in [p3|. 



1. Double counts 



2. Decay of propagators 



The integrals in ( |5^ are controlled by two mecha- 
nisms, the de-phasing given by momentum integrals and 
the decay of propagators. Let us assume for a while 



The double count of correlations is a serious mistake for 
theories which aim to go beyond the Boltzmann equation. 
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Here we show that the decay of propagators during the 
cohision removes this double count. 

We use the simplest approximation of the decay, 



-le 



(59) 



where = e{k,t - §) and 7^. = ^^(^k,t-^){k,t - §) is 
the pole value of 7 = ct-* + cr^ = —2lma^. The time 
argument shows that external or mean fields can mod- 
ify the propagation during the collision. The collected 
propagator thus reads 



= S„ 



^ _ _^ T- ) g-(7fc+7p+7fc-<j+7p+g)i 



(60) 



The kinetic equation holds only if the collision duration 
is short on the scale of the lifetime, (27)^^. In this limit 
we can expand the exponential decay, 

e-(7.+7p+7.-,+7p+,)5 = 1 _ (^^ + + + 7p+g)^- 

(61) 

Since we are dealing with the correction terms, we can use 
a simple approximation of the kinetic equation, dtPk = 
'~ IkP, to rearrange the contribution of 7's to its in- 
tegrand, 

T T c) T 

- l^{lk+lp + lk~q + Ip+q) = ^g^D - -Dap. (62) 

The last term, 



-Dsp = ^l{'^~Pp)Pk-qPp+q + {\~ Pk)(Tp Pk-qPp+q 
+ Pk){l~ Pp)(T^_ Pp+q + {l-pk){l~Pp)pk-qO. 



p+q 



-f^k PpC^^ Pk-q)^^- Pp+q) - PfcO-p (l-Pfe-g)(l-Pp+g) 



-PkPp(yl^„{i-Pp+q) - PkPpil-Pk-q)<y. 



P+9' 



(63) 



can be neglected leading to the three-particle contribu- 
tion to the scattering integral. Please note that such 
neglects have not been necessary in our suggested con- 
struction p[f] outlined in Sec. II. 

Within the linear approximation in time gradients one 
finds that the decay reduces the retardation, 

e-(7»=+7p+7fc-,+7p+,)i _r) ^ D (^t-^,Ty (64) 

Now the collective propagator and distribution have the 
same time retardation of the center-of-mass time being 
equal to ^ . This time shift removes the double count and 
reproduces the result obtained already in Sec. II. 

We note that the original Levinson equation |^2| does 
not include the decay of propagators in the scattering 
integral. The presented discussion shows that the de- 
cay plays the important role for the correct treatment of 
renormalizations and other contributions of correlations. 



B. Collision delay 

We have seen that within the Born approximation the 
retarded collision integral of the Levinson equation can 
be transformed into the instantaneous scattering integral 
of the Boltzmann type and quasiparticle renormalizations 
of the single-particle distribution and energy. Beyond 
the Born approximation, there is the additional retar- 
dation expressed by times tr^a in the selfenergy (|4^). 
This retardation is responsible for the virial correction, 
in particular for the correlated density ric- Since the re- 
tardation times Tfi,^A do not enter the final states of the 
collision, they cannot be removed in the above manner. 
On the other hand, the time integrals over tr^a can be 
eliminated so that the collision duration is expressed by 
an effective time-non-locality of the scattering integral. 
For simplicity of presentation in this section we assume 
that T-matrices and selfenergies are independent of the 
center-of-mass time. This will lead us to the basic idea. 
The complete result with all time and space dependences 
are given in the next Sec. IV. 

Within the linear approximation in time gradients we 
can treat the contributions of times tj^^^a separately from 
other retardations. We first implement the quasiparticle 
ansatz and expand in retardations. 



[ k ~ q,t - ^{tr + TA + t),t - TR + ta 
X g< [ p + q,t - 1-{tr + TA + t),t - tr, + ta 



/ 1 1 a 1 , ,9 



X fk-qit).fp+q{t)e- 



(65) 



The contribution of has been discussed above, now 
we focus on the term due to ^{tr + ta)^- 

Introducing the energy representation of the T-matrix, 



T^iTR) = J '^t-M 



(66) 



and a complex conjugate for T , the contribution of 
d_ 

2\'n I 'A/Qt 



]iTR + '''a)-§t to the selfenergy ([48[) reads 



(27r)e 

dijjRd(jJA,R,, .,Af-i N 

—pr-T^t [k,p,q,ujR)t (k,p,q,uJAj 



X J dTRclTAe-'^^^'^e"^^^^ {tr + ta) 

X ^/fe-«W-'(^'=-+^-+')("-"«+"^). (67) 

The integral over tr^a results in the derived (5-functions 
which can be readily integrated out. 



9 



A(T<(fc,t,r) 



dpdq 
(2^' 



2i \ duj duj 
d_ 



UJ — Sk-q+Sp+q 



(68) 



The correction (^8|) is proportional to the collision de- 
lay defined from the phase shift, 



duj 



(69) 



Indeed, writing the T-matriccs in the form t^'"^ — |t|e^*"^ 
where |ip is the amplitude of the scattering rate and ( 
is the scattering phase shift, we obtain 



2i \ duj duj 



(70) 



This allows us to join the correction with the time-local 
part of the selfenergy, 

dpdq 



a<{k,t,T) = 



(71) 



(27r)6 

X fk-g{t^ At)fp+g{t~ At) 



The selfenergy (|7l| ) can be compared with (|48|). First, 
the retardation by — r/2 has been suppressed because it 
is eliminated by the method described in Sec. II. Second, 
the retardations by the internal times of the T-matrices, 
tr^a, have been reduced from the integral form of ( ^ ) 
to the effective time non-locality of (|7^) . The latter form 
includes only a single characteristic time which is sim- 
ilar (but not fully identical) to the collision delay defined 
by Wigner from the scattering phase shift. Note that in 
form (|7^) no time integration appears. 

An interpretation of the collision delay At is troubled 
by the fact that this effective time does not have to be 
positive. Apparently, one cannot interpret the collision 
delay as a mean time scale of the tu^a integrals in the 
selfenergy (|4^). The negative value of At appears when 
the phase shift decreases with increasing energy, as it 
happens above a bound state or a resonant scattering 
state. For instance, the collision delay of nucleon-nucleon 
collisions is negative at low relative energies due to the 
deuteron bound state, see . The negative value of the 
collision delay merely tells that particles anti-correlate at 
a given relative energy. 



IV. NONLOCAL SCATTERING INTEGRAL 

The selfenergy ( |7l| ) is easily converted into the scatter- 
ing integral. The |tp represents the scattering amplitude 
in the T-matrix approximation. The exponential func- 
tion turns into the energy conserving (5-function. Com- 
pared to the common scattering integral of the Landau- 
Silin equation, the only difference consists in the collision 



delay At entering the initial states of the collision. Of 
course, in the previous section we have assumed a homo- 
geneous system with time-independent selfenergies and 
T-matrices. For a general case more nonlocal corrections 
appear. These corrections we discuss in this section. 

In this section we thus focus on the scattering inte- 
gral. Our aim is to derive the scattering integral within 
the same approximations as the skeleton kinetic equa- 
tion discussed in Sec. II, i.e., keeping all hnear gradients. 
This procedure leads to the non-instantaneous and non- 
local scattering integral. We use the Galitskii-Feynman 
T-matrix approximation of the selfenergy which is ap- 
propriate for dense systems of particles interacting via 
short range potentials. The T-matrix approximation for 
non-equilibrium Green functions is briefly presented in 
Appendix 0. 



A. Quasi-classical limit 

The consistent treatment of gradient contributions is 
provided by the quasi-classical limit, i.e., the linear ex- 
pansion in gradients. The collision delay derived above 
shows that beside gradient contributions forming the 
drift in the skeleton equation, there are non-trivial gra- 
dient contributions making the scattering integral non- 
local in time. Transforming the system into a running 
coordinate framework, one can see that the collision de- 
lay has be accompanied by the space non-locality of the 
collision integral. Accordingly, the gradient expansion 
has to be applied also to all internal space and time in- 
tegrations of the scattering integral. 

The quasi-classical limit and the limit of small scat- 
tering rates explicitly determine how to evaluate the 
scattering integral from the selfenergy a^. For non- 
degenerate systems, a very similar scheme was carried 
out by Barwinkel One can see in Barwinkel's pa- 

pers, that the scattering integral is troubled by a large 
set of gradient corrections. This formal complexity seems 
to be the main reason why most authors either neglect 
gradient corrections at all |]55| , |66| o r provide them buried 
in multi-dimensional integrals |30| , |3l| , |67| . For a degener- 
ate system, the set of gradient corrections to the scatter- 
ing integral is even larger than for rare gases studied by 
Barwinkel, see ^l[| . To avoid manipulations with long 
and obscure formulas, the gradient corrections have to 
be sorted and expressed in a comprehensive form. To 
this end, we will express all gradient corrections in the 
form of shifted arguments, as we have done above for the 
collision delay. 

The quasi-classical limit of the scattering integral with 
all linear gradients kept is a tedious but straightforward 
algebraic exercise, see Appendix ^ giving 



■(w,fc,r, t) 



10 



■ dp dE dq dn laAa 
(27r)3 2^ (27r)3 2^ I " 2~d^ 



xg>iE,p,r- A2,t) 



xg<{uj 



Q - Ae, k 



xg<{E + n- Ae,P' 



q ^ 
q 



AK,r ^ 



A3,t^ 

- A4,t- 



At) 
-A,). 



(72) 



With g^'^ substituted from the quasiparticle approxima- 
tion, this formula turns into the functional of the quasi- 
particle distribution. As one can see, the initial states, 
k — q and p + q, include the collision delay A(, as in the 
homogeneous system describe by (|7l|). 



Beside the collision delay, non-local selfenergy ( 72 ) in- 
cludes a number of other gradient corrections which enter 
nearly all arguments. These are expressed via A's which 
are derivatives of the scattering phase shift where 



R,A 



lile^F^-P, according to the following list |5l||5g] 

A2 
A3 

A4 



At = 



dp dq dk ' 

^dk' 
d(j) d(p 
dq dk ' 

duj' 
2 

~2 at ' 

] (A2 -t- A3 -t- A4 



(73) 
(74) 
(75) 
(76) 
(77) 
(78) 
(79) 



Finally, there is the energy gain which is discussed in [|68| . 
Numerical values of the delay Af and the displacements 
A2_4 for the scattering of isolated particles are presented 
in(§]. 

Sending all A's to zero one recovers the instantaneous 
and local approximation of the selfenergy, Eq. ( [z^ ) . The 
non-instantaneous and non-local corrections given by A's 
appear in three ways: in arguments of the correlation 
functions, in arguments of the T-matrix, and as a fore- 
factor. In correlation functions, the A's describe dis- 
placements of the initial and final positions of colliding 
particles, the final state of one of the particles is fixed to 
the balanced phase-space point {k,r,t). The arguments 
of the T-matrix correspond to the center of arguments of 
all initial and final states. 

The pre-factor has no special physical meaning as it 
depends on the actual choice of energy and momentum 
variables. For instance, using in ( [72[ ) an integral over 
E' = E — Ae instead of the integral over E, a corre- 



sponding factor, dE /dE' 
pears. 



1 + d^AE = 1 - ^dtAu ap- 



B. Extended quasiparticle approximation 

Now we can complete the derivation of the kinetic 
equation. Since the scattering integral is already pro- 
portional to the 'small' scattering rate, we can use from 
the extended quasiparticle approximation (|2^) the pole 
part, g^{uj, k, r, t) = /(fc, r, i)27r^(a; — e(fc, r, t)), to con- 
vert selfenergies a^'^ into functionals of the quasiparticle 
distribution. The selfenergy cr< is given by ( ^2[ ) and cr> 
is obtained by the interchange >< — ><. The resulting 
kinetic equation reads 

dfi _^ de^dfi _ dsi dfi 
dt dk dr dr dk 

= / ^^^'"^ -f £2 - £3 - £4 - 2Ae) 
1 dA2 de2 dA2 \ 



X 1 



2 dr dr dcu I 
t(ei+e2-AE,k-h\K,P~ ^Ak ,q,r-Ar,t~ ^A^ 



X ((1 - /i)(l - /2)/3/4 - /l/2(l - /3)(1 - /4)) 



The abbreviated notation of /'s and e's means 

ei = e(fc,r, i), 

£2 = e{p,r - A2,t), 

£3 = £(fc-q- Ax,r- A3,t- At), 

£4 = e{p + q - AK,r - A4,t - At). 



(80) 



(81) 



The bar reminds that all A's appear in arguments with 
negative signs. In the A's and their derivatives the energy 
uj + E = El + 62 IS substituted after all derivatives are 
taken. Note that the pre-factor has changed compared 
to ([72[). In (72), £2 depends on A2 which depends on the 
energy E. This E'-dependence has resulted in a norm of 
the singularity. 

The non-local kinetic equation ( |80| ) represents the 
main result of this paper. In the following we will discuss 
some symmetry properties of this equation. 



C. Markovian approximation 

In the scattering integrals of (^ij) the collision delay 
At appears together with displacements A2_4. Actual 
values of these types of corrections are closely linked. In 
the original Wigner's approach, the collision delay was 
identified from the position of the center of wave packets 
in the final state. Apparently, this position can be de- 
scribed either as the space displacement or via the delay. 
A corresponding rearrangement of the scattering integral 
is achieved using the free-space time evolution of the dis- 
tribution, e.g.. 
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h = f{k-q-AK..r-A3,t-At) 

~f{k-q- Ak + FsAt, r - A3 + fsAt, t), (82) 

where V3 = dk£z is a velocity and F3 = —drS^ is a force 
acting on the particle in the final state. Substituting the 
right hand side (and similar for /4 and £3,4) into (^l|), one 
obtains instantaneous scattering integrals. This liberty of 
rearrangement is only approximative, some quantities are 
sensitive to such approximations. For instance, the cor- 
related density discussed below is directly proportional 
to At as it measures an effective number of particles in 
the collision process. On the other hand, the rearrange- 
ment into the instantaneous form makes the kinetic equa- 
tion Markovian what significantly simplifies its numerical 
treatment, in particular if negative collision delay takes 
place. 



h = G>E< - G<S> and Ir = 5]<G> - E>G<. For the 
Galitskii-Feynman selfenergy, 

/; = G>G>T"G<G^T'^ - G<G<T"G>G^T'^, (84) 

we have to add and subtract GfGfT^G^GfT^ to the 
second term of (|^ to achieve the Galitskii-Feynman 
two-particle spectral function, 

/; = Ai2T'^G<G<T^ - G<G<T'^A34T^. (85) 



In the last term we can apply the optical theorem (|8^ 
for the interchange T^^34T^ — > T^A34T^, therefore 



/; = Ai2T'^G<G<T^ - G<G<T^A34T". 



(86) 



D. Space-time symmetry of non-local corrections °^ phase shift (j), 



In the same way one can rearrange Ir- 

Technically, the interchange of the retarded and the 
advanced T-matrices is equivalent to the change of sign 



t 



R _ 



Itle-*"^ and 



The scattering-out integral in (jsy) has the same sign 
of A's as the scattering- in. Accordingly, due to the non- 
local and non-instantaneous corrections, the scattering- 
out cannot be obtained from the scattering-in by the 
space-time symmetry, as it is common in the classical 
approach to the kinetic equation. This problem has a 
practical impact. The simulations of pressure with non- 
instant corrections has been discussed by j6^] . In Monte- 
Carlo simulations based on semi-classical trajectories the 
space-time symmetry is an unavoidable part of the scat- 
tering integral. To understand internal times of collisions 
and non-local corrections, it is thus necessary to clarify 
the space-time symmetry. 

It is possible to show that the scattering-in and - 
out are connected by the space-time symmetry, one has 
to take into account, however, the non-local and non- 
instantaneous corrections to the T-matrix itself. To this 
end we employ the optical theorem which can be ex- 
pressed in two forms. 



ImT = r^^r"^ = T^AT'^. 



(83) 



Here, symbol Im denotes the anti-Hermitian part. The 
multiplication includes integrals over two-particle states 
and times. The function ^34 = G3 G4 — G^G4 = 
G> (3, 3')G> (4, 4') -G< (3, 3')G< (4, 4') is the two-particle 
spectral function which includes the Pauli blocking of 
internal states in the Galitskii-Feynman approximation. 
Both forms of ImT can be derived by algebraic manipu- 
lations with the ladder equation for T^'^ and hold out 
of equilibrium. 

Implementation of the optical theorem (|8^) requires 
to rearrange the scattering integrals. This is most conve- 
niently done on the level of the starting non-equilibrium 
Green functions. The integrand in the right hand side 
of (^ is composed of two terms, in a reduced notation 



t'^ — \t\e^'^ . The implementation of the optical theorem 
thus has no effect on the local parts of the scattering 
integral which depend exclusively on the amplitude \t\ 
of the T-matrix. In contrast, all non-local corrections 
given by A's ( [73| - |7S| ) are linear functions of the phase shift 
<j) and thus reverse their signs. In the kinetic equation 
(^), the combination corresponding to ^^34 is achieved 
by regrouping distributions as 

(1 - /l)(l - h)hh - hh(l - /3)(1 - h) 

= (1 - /i - l2)hh - hh{i -h- h). (87) 

The intermediate-state Pauli-blocking factor 1-/3-/4 is 
the occupation number corresponding to ^134. Inverting 
signs of the phase shift, i.e., of A's, related to the scat- 
tering integral with the Pauli blocking 1 — /a — /4 one 
arrives at the kinetic equation 



dfi ^ dei dfi dei dfi 
dt dk dr dr dk 



dp dq 
(2^(2^ 



1 



1 9A2 de2 dA2 



2 dr 



dr duj 



1 



t{ei+e2-AE,k-]^AK,p-]^AK,q,r-Ar,t-^At 

X27:d (ei +82 -63 -Si- 2Ae) (1 - /l - /2)/3/4 

dp dq f 1^^2 de2dA2\ 
(27r)3 (27r)3 ^ 2 dr dr doj J 

t(ei+e2 + AE,k+^AK,p+^AK,q,r+Ar,t+^At 
x2T:d (ei +82-63-64 + 2Ae) hh{l ~ fs - h)- 



In functions with bars arguments are given by (^), in 
functions without bars arguments are identical except for 
the reversed (i.e., positive) signs of all A's. 
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The kinetic equation with the space-time symmetry 
of scattering integrals reminds to the Enskog equation 
for classical hard spheres with respect to two features. 



First, if one suppresses the collision delay, see (|82[), 
the scattering-out becomes the space inversion of the 
scattering-in. Second, the scattering rate (given by |tp) 
is centered between initial and final states of collisions. 

The kinetic equation (|8|) has a difficult interpreta- 
tion because its scattering integrals have positive or 
negative values following the sign of the Pauli block- 
ing factors 1-/1-/2 = (1 - /i)(l - /2) - /1/2 and 
1 - /3 - /4 = (1 - /3)(1 - h) - hh- This strange be- 
havior appears due to the stimulated transition processes 
given by the terms — /1/2 and —fsfi on top of the stan- 
dard Paufi blocking, (1 - /i)(l - /a) and (1 - /3)(1 - U). 
While in the local approximation of the scattering inte- 
grals the -in and -out stimulated transitions exactly com- 
pensate giving no net contribution, the non-local correc- 
tions translate this feature known from the equilibrium 
Green functions on the level of kinetic equation. Our 
difficulties with the classical interpretation of (|8|) thus 
may follow from our low understanding of stimulated pro- 
cesses in general. At very low temperatures, the stimu- 
lated transitions become so strong that the above quasi- 
particle picture breaks down and the system develops 
super-conductivity having an additional degree of free- 
dom and a different spectrum of single-particle excita- 
tions. Although the Galitskii-Feynman approximation is 
not suited to describe the super-conducting state, it re- 
veals a singularity at Tc which signals its onset. 

Let us put aside problems with stimulated processes 
and return to our discussion of characteristic times. For 
the moment we will assume that the collision delay At 
is positive. In the scattering-in integral which is the first 
term on the right hand side of (^8|), the initial condition 
is at the time t — Af, i.e., before the time t at which 
we balance the changes at the phase-space point {k,r). 
This has a clear interpretation, at t — At the particle of 
momentum k — q enters into the collision process from 
which it is released at t into {k,r). In the scattering- 
out process, the particle leaves the cell (fc,r) at t when 
it enters the collision process. This process will end at 
t + At when we have to check the accessible phase space. 
While the scattering-out has a natural interpretation, its 
description violates the causality. Indeed, the occupa- 
tion factors in future time t + At are supposed to de- 
cide whether the process comes through or whether it 
is blocked. We note that the violation of the causality 
on the microscopic time scale is common in kinetic equa- 
tions. Already the Fermi Golden Rule includes a time in- 
tegral into future which eliminates off-shell processes. In 
fact, it is the Fermi-Golden- Rule- type approximation, in 
this paper represented by the pole approximation, which 
makes negative values of the collision delay possible. 

For a detailed discussion of the space-time versus 
particle-hole symmetry for various forms of the two- 



particle spectral function see |7C 



E. Correlated density 

As mentioned in the introduction, in the system de- 
scribed on the level of the Galitskii-Feynman approxima- 
tion the density of quasiparticles, nf = f j0^f , differs 
from the density of composing particles n ^ J j0j^p by 
the amount called the correlated density ric = n — Uf. 
After some algebra, from relation (p2[) one finds 



dkdpdq ^ 



X(5(£i + £2 - £3 - e4)/l/2(l - /3 - fi)- 



(89) 



The functions |t|, e's and /'s are in the local approxima- 
tion having no A's in arguments. This expression is the 
generalized Beth-Uhlenbeck formula ||,|6|j4|j7|-^,|9| . 
The correlated density is proportional to the energy- 
derivative of the phase shift expressed here via the colli- 
sion delay. 

This correlated density is consistent with the equation 
of continuity found from the kinetic equation (p8|). For 
simplicity we will assume a homogeneous system for the 
equation of continuity reading, dtn = 0. Taking the mo- 
mentum integral over (88) one finds 



d_ 

dt 



dk 
(2^ 

X(5(£i + 



h = - 



d f dkdpdq 



t'At 



£2 



dt J {2tt)» 

£3 - £4)/l./2(l - /3 - /4) 



(90) 



Using ( |89| ) in the right hand side, one finds dt{nf +nc) = 
0. It should be noted that the number of quasiparticles is 
not conserved, dtUf ^ 0, because during the time interval 
At the colliding particles are excluded from the single- 
particle statistics being in two-particle scattering states. 

One can see that instantaneous approximations of the 
kinetic equation cannot capture the correlated density. 
In the same time, the correct description of the quasi- 
particle density makes the choice (|7^) of the collision de- 
lay preferable to other characteristic times one can define 
in the spirit of Wigner from the asymptotic behavior of 
final states of binary collisions. In Appendix ^ we com- 
ment on some mistakes done in previous studies of the 
correlated density. 

The complete proof of conservation laws with explicit 
expressions for the correlated pressure and energies can 
be found in ImI. 



V. SUMMARY 

In absence of the time operator, there are many defi- 
nitions of characteristic times one can associate with col- 
lisions. Regardless of a definition one uses, if the finite 
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duration of collisions is included in a kinetic equation, 
this equation is non-Markovian since the initial and final 
states of the collision are at distinct times. The family 
of non-Markovian kinetic equation also includes retarded 
equations of Levinson type in which the collision is ex- 
pressed in terms of the time integral describing the whole 
process of the two-particle interaction. 

For slowly varying systems, we have shown that 
one half of the Levinson-type retardation describes the 
off-shell motion and corresponding renormalizations of 
single-particle functions. The second half compensates 
the decay of propagators during the integration and can 
be eliminated leaving the pole contribution to the scatter- 
ing integral identical to the one obtained from the quasi- 
particle approximation. 

Finally, we have discussed the collision delay result- 
ing from the energy dependence of the scattering phase 
shift. We have shown that this non-Markovian correction 
is connected with other non-local corrections to the scat- 
tering integral which gives us a freedom to define the col- 
lision duration in different ways. The physical quantity 
sensitive to the actual choice of the collision delay is the 
density of quasiparticles, or its complementary quantity, 
the density of correlated particles. The method presented 
in this paper, and at the same time the choice of the col- 
lision delay, reproduces the correlated density obtained 
within the generalized Beth-Uhlenbeck approach. 

The non-instant and non-local corrections given by 
the A's do not change the structure and the overall in- 
terpretation of the scattering integral but only slightly 
renormalize its ingredients. The exclusive dependence of 
the non-local and non-instant corrections on the scatter- 
ing phase shift confirms results from the theory of gases 
29|J32 3^ obtained by very different technical tools. 



The non-instant and non-local scattering integral in the 
form ( ^sl ) parallels the classical Enskog's equation, there- 
fore it can be treated with numerical tools developed for 
the theory of classical gases, see e.g. Q at least as long 
as 1 — /i — /2 remains positive. 

The virial corrections to the balance equations appear 
from intrinsic gradients instead of correlation parts in 
the equation of the reduced density matrix. In |]52| it is 
shown that the nonlocal kinetic equation leads to com- 
plete balance equations for the density, energy and stress 
tensor which establish conservation laws including corre- 
lated parts. The nonlocal kinetic equation derived here 
requires no more numerical efforts than solving the Boltz- 
mann equation. The numerical solution and application 
has been demonstrated in l7a-f77 
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APPENDIX A: T-MATRIX APPROXIMATION 

In order to describe short-ranged two-particle interac- 
tions, it is necessary to introduce the standard approxi- 
mation of the many-body theory, the T-matrix approxi- 
mation In the following we give a brief com- 
pilation of important formulas. 

The sum of ladder diagrams is defined as the causal 
T-matrix 

(12|r|l'2') = V(12l'2')+V(1233')G(31)G(3'2)(12|r|l'2'). 

(Al) 

Arguments not present on the left hand side are in- 
tegrated over. For simplicity, we do not assume the 
sum over spin and isospin and the anti-symmetrization 
for identical particles. Since the potential is local in 
time, we can simplify the T-matrix as (12|T|1'2') = 
{xiX2ti\T\x'ix'2t'i)5{ti — t2)S{t[ — t'2), and ladder sum- 
mation (Al) reads 



{xiX2t\T\x'^x'2t') = V{xiX2x'-^x'2)6{t - t' ) 

+ V{xiX2Xzx'^){xzx'^t\Q\xiX2l){xiX2i\'T\x'ix'2t')^ (A2) 

where we have introduced the two-particle Green func- 
tion, 

{xiX2t\Q\xiX2'£) — G{xitXiT)G{x2tX2t). (A3) 

The causal selfenergy then reads 

I](ll') ~ ~^ j dX2dx'2G{x'2t'iX2ti){xiX2tl\T\x'iX2t'i). 

(A4) 

Using the Langreth-Wilkins rules (which arc equiv- 
alent to the algebra on the Kcldysh contour) we obtain 
the real-time Green's functions. From (A4) follows 



S^/-^ = -iQ{ti - t2)[T>G< + r<G>] 



r^G< 



T<G^. 



(A5) 



We have abbreviated the notation, all integrations and 
variables are the same as in (|A4|). From ( [A^ ) follows 
|5|,|6|j80) 



jRjA ^ -^gRIArj-RIA ^ 



(A6) 

(A7) 



Explicitly (A6) reads 
< xxX2t\T^\x\x'2t' > 

= J dxidx\dx2dx2dtdt' < xiX2t\%^\xiX2t > 

X < XiX2t\g^i,\x\x'2t' >< x\x2t'\T^i\x[x'2t' > 



(A8) 
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In the scattering integral we use the following repre- 
sentation of the two-particle T-matrix, 

/ ' W\ f dkdpdqdVt 

{xiX2t\T\x^X2t ) = J — (27r)To— ^ (^^ 9' ^) 

X exp [i {kxi + px2 — (fc — q)x[ — (p + q)x'2 — fl{t — t'))] , 

(A9) 

where r — j{xi + X2 + x[ + x'2) and t = ^{t + t') are 
center-of-the-niass coordinate and time. 



1. Quasiparticle approximation 

We introduce the notation for the known quasiparti- 
cle approximation without gradient corrections. In the 
quasiparticle approximation the free two-particle propa- 
gator Q (denoted in the mixed Wigner representation as 
C) takes the form 

The selfenergy then reads 



CT^ (uj, k, r, t) 



t>{u; + e{p,r,t),k,p,0)f{p,r,t), 



(27r 



(All) 



and analogously a^. In the rest of this chapter we sup- 
press variables r,t. From (A7) one finds 



t<in,k,p,o) 



(27r)3 



\t <| {ek-q+ep+q,k,p,q) 



X fk-qfp+q2n6{uj - Ep+q ~ £fc_g), (A12) 

and is given by interchange / ^ 1 — /. 

For the determination of the quasiparticle energy one 
needs the retarded selfenergy 



tT-"(w,fc) 



dp duj' (lu' + Ep, k,p,0)fp 
(27r)3 2tt lo' — bj — irj 
dp diu' t<{Lu' + ep,k,p,0){l- fp) 



-I 



(27r)3 27r 
dp 



Uj' ^ UJ — irj 
t^iu + ep,k,p,0)fp 



(27r)3 

dpdq \t^{ek-q + ep+q,k,p,q)\'^fk-qfp+q 



(27r)6 



£k-q + £p+q - Sp - U + IT] 



(A13) 



In the real part of the selfenergy the real part of the 
T-matrix appears, 



Retain, k,p,0) 



dq 1 — fk-q ~ fp+q 
(27r)3 n ~ Hk-q ~ f-p+q 

X \t^\'^{ek-q + ep+g, k,p, q). (A14) 



This form follows from the Kramers-Kronig transforma- 
tion oit> - t<. 

The ladder equation without gradients reads [notation 
is specified by (B2)] 



{Pl\t'\q')\p2) ^ {pi\V{q')\p2) + 



Up'dp" 
(2^ 



y.{pi\V{q')\p'){p'\C^{q')\p"){p"\tfMp2) 



(A15) 



In the coordinates, k = ^ + Pi:P — \ ^ PIt 1 ^ Pi ^ P2, 



corresponding to (A9) this equation takes the form 
t'^{n,k,p,q)^V{k,p,q) 



dp'dp" 



v{k,p,p'K^{n,k-p\p + p\p' -p") 



(27r)6 

xt«(f],fc-p",p + p",g+p") 



(A16) 



Since ^^(fc,p, g) = (27r)3(5(q)(^^(fc)(7^(p), in the quasipar- 
ticle approximation one obtains 



t^{n,k,p,q)^V{k,p,q) + 



dq' 
(2^ 



V{k,p,q') 



1 fk-q' fp+q' j^Rlm I ' I ' I '\ 

X — . t ^{il,k - q ,p + q ,q + q ). 



^ - £p+q' - £k-q' + iri 



(A17) 



APPENDIX B: GRADIENT EXPANSION 

Let us consider a product of two-particle functions 

C(1234) = J d3'd4'y^(123'4')6(3'4'34). (Bl) 

We transform this product into mixed representation 
( |A9| ) keeping gradients till the linear order. 

Variables in ( [A9| ) follow variables of the scattering in- 
tegral. For the scattering-in, the momenta correspond to 
the process k — q,p + q — > k,p. These variables, how- 
ever, do not form a convenient algebra for the gradient 
expansion. To this end we use the coordinates 

a = 1 - 2 

/3 = 3-4 



r= 1(1 + 2-3-4) 



x== -(1 + 2 + 3 + 4) EE (r,t). 



(B2) 



Since the relative coordinates, a and /3, obey the stan- 
dard matrix algebra while only the coordinate r and x 
undergo the gradient expansion, we write these argu- 
ments in form {a\C{T, x)\P). In this notation the product 
(Bl) reads 
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{a\C{T,x)\l3) = / d-idf{a\A ( r - f, a; + -f ) I7) 



X {l\B ( T,x+ -(r -r) 



(B3) 



Assuming a slow variation of the center-of-mass coordi- 
nate X we can write this product keeping gradients in x 
up to the linear order. The relative coordinates can be 
suppressed being independent of this expansion. Via the 
Fourier transformation of coordinate r, 



C{t,x) 



dn 
(2^ 



(B4) 



we express the operator C as a function of the sum mo- 
mentum, K2, K3) — k + p, and energy K4 = O of both 
particles. One can see that the product has the same 
form as products of single-particle functions, 



2 I % 
c = ab [a,b] = ah I 1 [In a, In 6] 



(B5) 



where the (generalized) Poisson bracket applies to the 
sum coordinate and time. The logarithmic form of (B5) 
is merely a convenient notation because a and b are ma- 
trices in a and p. The logarithm is taken of each matrix 
element of a and b and summed together with the prod- 
uct ab. 

In parallel with the Dyson equation, there are no linear 
gradients in the ladder equation (A7). Like the retarded 
Green function, the T-matrix is given by the matrix in- 
version, Tj^^ = — Q^- Since the T-matrix is symmet- 
ric with respect to matrix arguments, {s\T\s') — (s'|T|s) 
(s and s' are momenta associated with a and /?, respec- 
tively), the matrix inversion does not bring any gradients. 

The selfenergy, S<(13) = T<(1234)G>(42), includes 
a number of gradient contributions due to internal gra- 
dients in and a simple non-local correction due the 
convolution with . 



1. Two-particle matrix products 



For the calculation of the selfenergy in (B15) we need 



t< ( r!,/c,p,0,r- iA2,t 



that all gradient corrections depend exclusively on the 
derivatives of phase shift (j)- Using definitions ( |73[]79| ) 
and the linear approximation, a{x){l — lS.dxa) = a(a;— A), 
expression ( p36| ) can be given a form 

dq dQ 
(2^(2^ 

xC<(r!-2A^,/c-g-Afe,p+g-Afc,Q,r-^-^,t-At 
xi^ (r!-A„,fc-^,p-^,g + Q,r-A„i-^ 



(B7) 



for the zero transferred momentum, (7 = 0, and its 
infinitesimal vicinity. Making the gradient expansion of ^ ^^^P 

[~i{k ~ q){xi - x[) - i{p + q){x2 ~ x'^)] 



(A6) in terms of the Poisson brackets one finds 



t< = iV<^( 1 - ^([lnt«,lng<] - [lnt^,ln5<] 



Factors of half in the non-local corrections to the ampli- 
tude result from the fact that this amplitude enters 
the collision integral in the square. 

In the absence of gradients, the two-particle function 
is singular in the 'transferred' momentum Q giving the 
only contribution for Q = 0. When the gradient correc- 
tions are already in the explicit form, one can employ this 
symmetry as it has been used above. Indeed, the com- 
plex conjugacy of the advanced and retarded T-matrices, 
= , requires equal arguments of both functions. 
Now we show that the infinitesimal vicinity of Q = 
brings additional gradient contributions. 

2. Convolution of initial states 

The two-particle correlation function representing 
initial states of the collision in (B7) is not a suitable in- 
put for the kinetic equation. We have to express in 
terms of the convolution of two single-particle functions. 
This convolution brings gradient corrections by which ef- 
fective positions of particles entering the collision process 
become distinct. 

^From the inverse transformation to ( [A9| ) we obtain 

C<{n,k-p,p + q,Q,r, t) 

= 2"^ / dxidx2dx\dx2dtdt' exp[iQ{x[ — xj,) + iil(t — t')] 



(B8) 



xS{xi +X2 + x[ +X2~ 4r)(5(t + t' - 2t) 
X G< {xit, x[t')G< {x2t, x'^t') . 



The substitution, xi = r + a + /3i/2, x'l — r + a — f3i/2, 
X2 = r - a + /32/2, x'2 = r - a - ^02/2, t = t + t/2, 
t' — t — f/2, shows that the (5-functions yield f = t and 
a — a and we obtain 



|ipg<(l -[(/), In .g<] -[(/), In (B6) C'^{n,k - p,p + q,Q,r,t) ^ 2^J dPidP2dadT exp[2iQa] 



In the second line we have decomposed the T-matrices 
into amplitudes and phase shifts, t^'^ = \t\e^^'^ . Note 



X exp [-i{k -q + Q/2)Pi -i{p + q- Q/2)/32 + iflr] 
xg<{(3i,r + a, t, T)g< {132, r - a, t, r), (B9) 
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where the representation ( [To| ) has been used. 

Now we hnearize the center-of-mass dependence of 's 
in a. The factor a can be represented by a derivative in 
front of the integral, a — — jdq. The remaining 
integration over a results in Tr^d{Q). When substituted 
into (B7), the differential operator in front of the integral 
is treated with the help of the integration by parts giving 



(BIO) 



and Q is integrated out. In other words the a in the 
arguments of g in ( B9 ) can be replaced by nonlocal shifts 
valid up to linear orders in gradients. This leads with 



< (n,k,p,0,r- 



xg<{k 
xg<{p- 



q, r 
q,r 



dq 

Afc Afc 
A3,t- At,rj-2A^) 



Ar,t- 



(Bll) 



3. Convolution of T- matrix and hole Green function 

Next we evaluate the convolution of the T-matrix with 
the hole Green function which is required for the selfen- 
ergy. Writing the selfenergy in the above representation, 
[t = 1 - 3 and a; = 5(1 + 3)] 



S<(T,a;) 



dadl3G> { -t - (3 + a,x - ^{a + (3) 



-{a- l3),x- -{a + (3) 



(B12) 



we see that the convolution couples matrix arguments a 
and P with the center-of-mass variables. By substitution, 
A = ■^{a + P) and ^ = a — f3, and expansion in gradients 
we obtain 



T,^{t,x)= j dfi G'' — t,x) 

dfi {}! — T, x) 



1 



T ,X 

2' 



A-^ 
2 



x/.AA(A + ||r<(.:|,.)|A-|) 



(B13) 



The second and third terms are gradient corrections due 
to the convolution. 



Now we can transform the selfenergy (B13) into the 
mixed representation. The momentum representation 
of the matrix algebra is introduced via unity operators, 
e.g., 1 = / with (A - /i/2|s) = e-(^-^/2) and 

(s'|A+/i/2) = e"**'(-^+'"/2). For the non-gradient part the 
integration over A results in the known fact that only the 
diagonal element, s = s' , contributes. For the gradient 
contribution, A —{i/2){ds — ds'), and the diagonal el- 
ement is taken after the derivatives are performed. Since 
we have considered the gradient corrections to the T- 
matrix already in the last chapter and since the shifts 
are additive, it is sufficient here to use the zero-order 
approximation of , 



(.|t<|.') 



dsds 



{sr\s){s\(:<\~s){s\t''\s'), (B14) 



where all functions have the center-of-mass argument 
(k, a;). From t^'^ — |f|e^*"^ and condition s = s' we find 
that A —ds(j)- By the substitution into variables of 
the kinetic equation, (g,0) = s — s', {k,Lu) — k/2 + s, 
(p, 0) = k/2 — s and {r,t) = x, one confirms that 
—dgcj) — 5A2, see ([73|), and the selfenergy reads 



(^^{k,r, t) 



f dpdn . , ^ , f 1 dA2\ 
J j^,9n-AP,r-A„t) [l - ) 



X t< ( n,k,p,0,r-lA2,t 



(B15) 



In this expression, abbreviates the right hand side 



of (B14), because the non-local correction A2 is defined 
only with respect to the integral over internal states of 
the collision. The norm term, 1 — ^drA2 results from the 
interchange, dr{A2A) = {drA2)A + A2drA, one has to 
make before the derivative is expressed in terms of the 
displacement, e.g., ff^(r) — A2drg^{r) = g^{r — A2). 



Together with ( Bll ) we obtain the result ([t^). 



APPENDIX C: COMMON MET PITFALL WITH 
THE CORRELATED DENSITY 

We note that one of us made a mistake deriving the 
correlated density directly from the retardation in the 
Levinson-type equation ||5^. The expression found in 
[ p9[ is identical to the correlated density except for two 
points, there is a wrong sign and the amplitude is as twice 
as large. Since this is a pitfall often met in literature we 
want to give some details here that it can be avoided in 
future. 

The simple expansion of the Levinson equation (^ up 
to the first order in memory ||59| yields 



d_ 

dt 



P - 



d_ 

dt 



n. 



(CI) 
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To obtain balance equations for the density including 
virial corrections, one has to integrate (CI) over momen- 
tum k, sec ||5^,^. Following the arguments developed 
for the Boltzmann-type equations, in | |59t , the wrong as- 
sumption was used that the correlated part of the colli- 
sion integral, fic — J dkR, combines with the left hand 
side of (pi]), nj — J dkp, to establish the density conser- 
vation, -^nj + he) = 0. The correlated density derived 
in this way reminds the result known from equilibrium 
[|6|,||j7l]-|7|l. Note that in this picture, the Wigner dis- 
tribution in left hand side of the Levinson equation is 
treated as the quasiparticle distribution what is a mis- 
take done usually in density operator studies js^. The 
wrong sign follows from this last misinterpretation. By 
definition the momentum integral over the Wigner distri- 
bution yields the full density, J dkp = n. The conserva- 
tion law thus tells that either J dkR = or the Levinson 
equation does not conserve the number of particles. 

As we have shown in this paper the Levinson equation 
contains additional gradient terms which exactly com- 
pensate the exphcit time gradients (used to derive fif). 
The above double count follows from the double count of 
the off-shell contributions discussed in this paper. 

Finally, the tr^a retardation was neglected. Since the 
correlated density (|8^) found from the off-shell contri- 
bution to the Wigner distribution is consistent with the 
conservation of the number of particles found from the 
Tfi^A retardation, one can see that these two contributions 
cancel and the collision integral of the Levinson equation 
conserves the number of particles, as it should. 
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